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Abstract 



We show that conformal maps of simply connected domains with an analytic 
boundary to a unit disk have an intimate relation to the dispersionless 2D Toda 
integrable hierarchy. The maps are determined by a particular solution to the 
hierarchy singled out by the conditions known as "string equations". The same 
hierarchy locally solves the 2D inverse potential problem, i.e., reconstruction of the 
domain out of a set of its harmonic moments. This is the same solution which is 
known to describe 2D gravity coupled to c = 1 matter. We also introduce a concept 
of the r-function for analytic curves. 
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1 Introduction 



By the Riemann mapping theorem, any simply connected domain on the complex plane 
with a boundary consisting of at least two points can be conformally mapped onto the 
unit disk. However, the theorem does not say how to construct the map. We argue that 
the map is explicitly given by a solution of a dispersionless integrable hierarchy which is 
a mult i- dimensional extension of the hierarchies of hydrodynamic type 0. 

In this paper, we only consider simply connected domains bounded by a simple ana- 
lytic curve. We show that in this case conformal maps are given by a particular solution 
of the dispersionless 2D Toda hierarchy (see 0, || and references therein). Surprisingly, 
this solution obeys (and is selected by) the string equation familiar in topological gravity 
and matrix models of 2D gravity 0, |^ . 

One may characterize a closed analytic curve by a set of harmonic moments of, say 
the exterior of the domain surrounded by the curve and its area: 

Ck = - 1 I z'^dxdy, k>l; Co = f f dxdy (1.1) 

J Jexterior J ^interior 

Here z = x + iy and it is implied that the point z = is inside the domain, whereas the 
point z = oo is outside. The integrals at A; = 1, 2 are assumed to be properly regularized. 
We assume this (infinite) set to be given and address two problems: 

(1) To find the harmonic moments of the interior 

= [ [ z^dxdy, > 1, (1.2) 
J ^interior 

and to reconstruct the form of the domain. This is the 2D inverse potential problem (see 
e.g. [§]). It is known that the solution, if exists, is unique if the domain is star-like. It is 
not known whether this is true for arbitrary domains with a smooth boundary. However, 
if the boundary is smooth, then a small change of moments uniquely determines a new 
domain (also with a smooth boundary). The integrable hierarchy gives the solution to 
this problem through the r-function, thus suggesting a new concept - the r-function of 
the curve. 

(2) To construct an invertible conformal map of the exterior of the unit circle onto 
the exterior of the domain. 

Explicitly, the area and the set of moments of the exterior (and their complex conju- 
gate) are identified with the times of the dispersionless 2D Toda hierarchy as follows: 

t = to = — , 4 = ^, 4 = ^, k>l. (1.3) 

These parameters should be treated as independent variables. We prove that the deriva- 
tives of the harmonic moments of the interior C-k with respect to the tj, ij enjoy the 
symmetry 

dC-k _ dC^j dCLk_dCLj_ 
dtj dtk ' dtj dtk 
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Therefore, the moments of the interior can be expressed as derivatives of a single real 
function which appears to be the logarithm of the r-function r(t, {tfc}, {tk}) of the dis- 
persionless 2D Toda hierarchy: 

dlogr _ C^k dlogT 
Vk = = ^7 — ' ^fc = = r.T ' k>l (1.4) 

71 Otk 71 Otk 

The equation which determines the r-function (the dispersionless Hirota equation) and 
therefore the moments of the interior is given in Sec. 5. Supplemented by a proper initial 
condition provided by the string equation, it solves the inverse moments problem for small 
deformations of a simply connected domain with analytic boundary. The r-function at 
tfc = 0, /c > 3 is given in the Appendix. 

In its turn, the conformal map is determined by the dispersionless limit of the Lax- 
Sato equations of the Toda hierarchy. Let z{w) be a univalent function that provides an 
invertible conformal map of the exterior of the unit circle |w| > 1 to the exterior of the 
domain. Let us represent it by a Laurent series 

oo 

z{w) = rw7 -|- ^ UjW~^ (1.5) 

j=0 

The conditions that w = cx) is mapped to ^ = oo and r is real fix the map, so the 
potentials r, uj are uniquely determined by the domain. The potentials of the conformal 
map obey an infinite set of differential equations with respect to the times (harmonic 
moments). They are evolution Lax-Sato equations of the hierarchy. The series z{w) is 
identified with the Lax function. The Lax-Sato equations are derived in Sec. 4. 

The coefficients of the conformal map have an equivalent description in terms of the 
r-function: the inverse map w{z) is explicitly given by the formula 




logM; = logz- -— + V— 7^-7— logr (1.6) 



Two comments are in order. Many equations of Sees. 3-6 can be found in the literature 
on dispersionless hierarchies [0, H, H, |^, 2D gravity 0, |[ and topological field theories 
||10| , \n\ without reference to conformal maps and inverse potential problem; we adopt 
them and their proofs for conformal maps. This work is the development of our previous 
joint work |12| with Mark Mineev-Weinstein, where some of the results described below 
were applied to the problems of 2D interface dynamics. 



2 Analytic curves and the Schwarz function 

Let a closed analytic curve be the boundary of a simply connected domain on the complex 
plane with coordinates z = x + iy,z = x — iy. The equation of the curve can be written 
in the form 

z = S{z) (2.1) 

Thus the curve determines the function S{z) which is defined also outside the curve. 
Moreover it is an analytic function within some domain containing the curve. We call 
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the function S{z) the Schwarz function of the curve (see e.g. [0). Its singularities 
encode the information about the curve. 

The Schwarz function may not be an arbitrary function of z. It obeys the unitarity 
condition: as is seen from (|2.1|) , its complex conjugate function^ is equal to the inverse 
function, i.e., 

SiSiz)) = z. (2.2) 

Being an analytic function in a strip containing the curve, the Schwarz function 
can be decomposed into a sum of two functions 5'(z) = S'^^^{z) + S'^~^[z\. one, 
is holomorphic in the interior of the domain, the second, S**^"-*, is holomorphic in the 
exterior. Their expansions around z = and z = oo determine the harmonic moments 
in the exterior ( |1 . ![ ) and interior ( |1.2[ ) respectively: 



7r 



fc-i 



k=l 



TT 



EC- 

k=0 



kZ 



-fc-1 



This follows from the contour integral representation of the moments 



(2.3) 



Ck = - 
2i 



z ^S{z)dz, 



curve 



-oo < A; < oo 



(2.4) 



obtained from ( |1.1| , |1.2| ) with the help of the Green formula. 

The Schwarz function is closely related to conformal maps. Consider the conformal 



map ( |1.5| ) of the exterior of the unit circle to the exterior of the domain. The inverse 
map sends a point z to a point w{z). Let us invert this point with respect to the circle, 
w — {w)~^ and map it back. (It is known that the conformal can be extended into 
some domain through the analytic curve. We assume that the inverted point belongs to 
it.) This operation is carried out by the conjugate Schwarz function S{z) = z{{w)~^) or 
S{z) = z{w~^). Obviously, if z belongs to the curve, then S{z) = z. Therefore, we can 
write S{z{w)) = z^w"^). 



To summarize, we have two formal Laurent series (see (|1.3| , |1.4| ) for the notation): 



Siz) = Y.kt 



kZ 



k=l 



- + ^VkZ 



-fc-1 



(2.5) 



fc=i 



S{z) 



^ ktkZ 

k=l 



fc-1 



-fc-1 



(2.6) 



fc=i 



They are connected by the unitarity condition (|2.2|) . The latter is resolved by the con- 
formal maps 

z{w) = S{z{w~^)) (2.7) 
z{w~^) = S{z{w)) (2.8) 
In their turn, the conformal maps are given by the series^ 



ziw) 



rw 



+ E 



UjW 



j=0 



(2.9) 



comment on the notation: given a Lament series f{z) — we set f{z) = fj^^- 

^To avoid confusion, let us stress that the functions z{w) and z{w~^) are complex conjugate only on 
the curve, i.e., at \w\ = 1. 



4 



z{w-'^) = rw"^ + J2 UjW^ (2-10) 
i=o 

They and the unitarity condition ( |2.7| ) estabhsh relations between harmonic moments 
of the exterior tk, harmonic moments of the interior Vk and the coefficients of the confor- 
mal map uj. Below we find an infinite set of differential equations, which determine the 
evolution of the potentials and the moments of the interior in times - moments of the 
interior. 



3 Symplectic structure of conformal maps and gen- 
erating function 

Deformations of the domain and, therefore, of the conformal map reveal the symplec- 
tic structure. In this section we show that the pairs logw, t and z{w,t), z{w~^,t) are 
canonical (see also [0): 

{z{w,t),z{w~\t)} = l (3.1) 

where {, } is the Poisson bracket with respect to w and the area t (all moments of the 
exterior tk, h are kept fixed). For any two functions f{w,t), g{w,t) the Poisson bracket 
is defined by 

We refer to ( pTT| ) as the string equation^. This equation suggests that deformations of 
conformal maps form a group with respect to the composition. To prove it, let us rewrite 
the l.h.s. in two equivalent forms with the help of ( |2.7| ), ( p.8| ). First, let us compute the 
derivatives of z^w'^jt) in (|3.1|) treating it as the composition S{z{w,t),t) as is in ( p. 81) : 

dz{w~\t) _ dS{z,t) dS{z,t) dz{w,t) dz{w-\t) _ dS{z,t) dz{w,t) 
dt dt dz dt ' dw dz dw 

As a result, we obtain: 

= (3.4) 

Similarly, treating z{w,t) in ( p.lj) as the composition S{z{w^^,t),t) we get 

{,M,,(^-',t)} = -JJ<^^. (3.5) 

In the r.h.s. of these equations the derivatives in t is taken at fixed z or z and then 
understood as functions of w. 



Now, using the series ( pl5| ) and (|2.9|), we conclude that the r.h.s. of ( |3l^ ) is 1 plus 
positive powers in w. However, the series ( p.6|) and (p.lO|) tell that the r.h.s. of (|3.5| ) is 



1 plus negative powers in w. This prompts us to (|0| ). 



^See Sec. 7 for the history of this equation and references. 
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The rest follows from the symplectic structure by treating deformations of the con- 
formal map (2^) along the lines of the multi-time Hamilton- Jacobi formalism [|], |14|. 



Let us introduce the generating function Q{z,t) of the canonical transformation 
{logw, t) — > {z, S). Its differential, dQ = Sdz + logwdt, implies that 

S = ^^. \ogw = ^^ (3.6) 
Using the Laurent series for the Schwarz function (|2.5|) , we get 

n{z, t) = Y: tkz' + t\ogz- -voit) - E ^z-', (3.7) 
fc=l ^ k=l ^ 

where Vq obeys 

dtvo = 2\ogr (3.8) 

The integration constant Vq can be interpreted as a logarithmic moment (see ( p.l5|) 
below) . 

As the Schwarz function, the generating function can be represented as the sum of 
functions, whose derivatives, S'^^^(z) = dz^'^^^z)^ are analytic in the exterior and the 
interior of the domain respectively: Vl{z) = il^^\z) + Q^~'^{z) —Vq/2. They have a simple 
electrostatic interpretation. Say, QS^\z) is the (complex) 2D Coulomb potential in the 
interior of the domain created by a homogeneously distributed charge in the exterior. 
In its turn, the QS~\z) is the complex Coulomb potential in the exterior created by a 
homogeneous charge in the interior: 

fiW(^) = -ll . log(l - -)dx'dy' = -f:%' (3.9) 
vr J J exterior ^ z' ^ k=i ^ 

fi(-)(^) =- \og(z-z')dx'dy' = ^\ogz--y^ ^z-' (3.10) 

TTJJinterior ^ ^ tt ^fc=i 

The real and imaginary parts of the generating function treated as a function of a 
point z on a curve have a clear geometric interpretation. The real part is half of the 
squared distance from the origin to the point 2 on a curve while | {Q{z, t) — Q{zi,t)) 
is the area A{z) (counted modulo nt) of the interior domain bounded by the ray ip = aigz 
and a reference ray ip = arg2;i. So, if z is on a curve, we can write 

Q{z) = ^Izl"^ + 2iA{z) , (3.11) 

where the right hand side is defined up to a purely imaginary constant. Indeed, integrat- 
ing by parts Q{z) = S{()d( = ^zS{z) + | J"^ {S {()d( — (dS {()) and taking into account 
that dA{z) = -^^{zdz — zdz) and S{z) = z on a curve, we obtain ( |3.11| ). 



It follows from ( p.ll| ) that if z belongs to a curve, then the variation of the VL[z) with 



respect to the real parameters. 



= dtnSt + [{dt,^ + dt,^) St^ + i{dt,n - dt^n) Im Stk] (3.12) 

k=l 
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is purely imaginary. The proof is simple. Let us analytically continue the equality 



(the real part of (|3.11 )) away from the curve, 



n{z) + n{S{z)) = zS{z) 



(3.13) 



and take the partial derivative with respect to tj at a fixed z. Note that the tj-dependence 
of the second term comes from both Q and S. We get: 



+ 



dm 



<:=siz) 



dSjz) 
~dt~ 



Applying the unitarity condition (|2.2|) , we observe that the second term in the left hand 
side is equal to the term in the right hand side, so they cancel. Restricting the equality 
to the curve again, we conclude that 



dn{z) ^ dQ{z) _ ^ 



dU 



(3.14) 



whence all terms of the infinite sum in the right hand side of ( |3.12| ) are imaginary. 
Similarly, for the derivative with respect to t (|3.14|) reads dt^^ + dtfl = (since t is real), 
so the first term in (p.l2|) is also imaginary. 

As is mentioned above, Q is defined up to a purely imaginary 2;-independent term 
(which may depend on the tj). However, its variation with respect to rea/ parameters, as 
in ( |3.12| ), is also purely imaginary. Taking into account ( |3.8| ), we fix it by the condition 
that Vq is real. From (|3.9|) , (|3.10|) and ( |3.11|) it easily follows that 



Vq = — I log \z\ dxdy 

vr Jinterior 



(3.15) 



(the logarithmic moment). To prove this formula, one should continue the equality 
Vq — 2TZe^l^^\z) = 2Tleil^~\z) — (which is just the real part of (|3.11| )) harmonically 
away from the boundary into the interior and evaluate both sides at z = 0. 



4 Lax-Sato equations for conformal maps 

Let us now vary higher harmonic moments 1^. The differential of the generating function 
changes accordingly [0: 

oo 

dn = Sdz + logwdt + ^{Hkdtk - Hkdh) (4.1) 

k=l 

where 

are Hamiltonians generating the higher flows in tj:,tk- Here the derivatives in are 
taken at a fixed z that belongs to the curve. The second equahty follows from (|3.14| ). 
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The Hamiltonians, being treated as functions of z and t obey the integrabihty condition 

(4.3) 



Computing the derivative in (f4.2|), we obtain: 



i7=z^--^-y^— (44) 

and similarly for Hk- 

The Hamiltonians Hj written in terms of the canonical variables w , t determine the 
evolution of the conformal map z{w) and z{w~^) with respect to the hierarchical times 
tk (harmonic moments). These are Lax-Sato equations: 

^ = W,.W} (4.5) 

^ = {ff,. (4.6) 

Note that these formulas can be extended to j = where Hq = logw generates the 
flow to = t. Consistency of these equations acquires the form of the "zero curvature" 
condition dt-Hi — dt^Hj + {Hi, Hj} = 0. It is equivalent to ( |4.3|) rewritten in the variable 
w. 

Now we are ready to prove that the Hamiltonians Hi and Hi have the form 

1 



H,iw) = [z\w))^ + -[z^iw))^ (4.7) 



Hj(w 



{w) = {z\w-'))_ + l{z^{w-'))^ (4. 



The symbols {f{w))± mean a truncated Laurent series, where only terms with positive 
(negative) powers of w are kept, {f{w))o is a constant part (w°) of the series. 

To this end, let us differentiate Hj by w and t, and express the result in terms of 
derivatives of z and S. We begin with the formula 



dHj dz{w) z{w ^) dz{w) z{w 



-1^ 



dw dw dtj dtj dw 



which is a simple consequence of the definition ( [4.21) . Replacing here z{w) by S{z) as 
before, we get: 

dH^ ^ mdsizi 

dw dw dtj 

Using ( |2.6| ), we find that dt^S{z) = Y^kLidtjVkZ~^~^ is a regular function of w at w = 
and its Taylor expansion starts from w'^. Thus Hj is also a regular function in w. 
Moreover, from the (O) we find that Hj is a polynomial in w of the degree j. Thus, 
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being so, we have Hj = {z^)+ + {z^)o — \dt Vo. To complete the proof, let us find the 
w^-teim in the Laurent series 

dHj _ dSjz) dzjw) dz{w) dSjz) 

~dr ~ ~di dtj di dt~ ^ ■ ' 

It comes from the first term of the r.h.s. of this expression and, together with ( |3.7| ), 
gives the desired result: 2{dtHj)o = 2dtjlogr = dtdt^vo = dt{z^)o. The dynamics of 
the conformal map with respect to tk can be obtained from ( [4.5| , [4.6|) by the comlex 
conjugatiion. Note that the Poisson bracket changes the sign as w ^ w = that is 
just consistent with the minus sign in ([4.2|) . Hence Hj = Hj{w~^) as is in ([4.8|) . 

The Lax-Sato equations ( [4 .51 , [4.6| ) with the Hamiltonians ( |4.7| , [4.8| ) imply that the co- 
efficients of the conformal map obey an infinite set of non-linear differential equations 
known as the dispersionless 2D Toda hierarchy. The first and the most familiar equa- 
tion of the hierarchy is the long wave limit of the Toda lattice equation: d"^ [r"^) / dt"^ = 
dHog{r^)/dtidh. 

We also mention other relations between the conformal map and the Hamiltonians: 

oo 

z{w) = Hi + -h + E kikHk-1 (4.11) 

^ k=2 

oo 

z{w)z{w-^) =t + 2neY,ktkHk (4.12) 

k=l 

They can be immediately obtained from the unitarity condition ( |2.7] , |2.8| ) by comparing 
the positive and constant parts of the Laurent series in w and using ([4.7|, ^l8|) . These for- 
mulas illustrate an important relation between harmonic moments and conformal maps. 
If all moments tfc = for k > N, then the Laurent series of the conformal map (|1.5| ) is 
truncated: uj = for j > N. This, in particular, gives another proof of Sakai's theorem 



| 15| : if all but the first three moments t, ti, t2 of the complement of a simply connected 
domain on the complex plane are zero, then the domain is an ellipse. 

The Lax-Sato equations (^75|,|4.6|) supplemented by the string equation ( p.l|) give 
the complete set of differential equations for the potentials Uj of the conformal map as 
functions of moments t,tk,ik- They uniquely determine small deformations of the map 
and eventually the curve. 



5 Symmetry of moments and the r-function of con- 
formal maps and curves 

The integrable hierarchy suggests a concept of r-function for curves and conformal maps. 
The r-function solves the problem of moments, i.e., restoration of the moments of the 
interior Vk out of the moments of the exterior {t^} and the area t. We define the r- 
function r(t; {ti}, {U}) as a real function, which determines the moments of the interior 
by the formulas 

d\ogT d\ogT dlogT 

^^ = ^7 ' ^fc = or ) ^0 = o , (5.1) 

otk otk at 
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The very existence of the r-function is due to the fundamental symmetry of harmonic 
moments: 

rill- rlu, rill- rlii, riDr, FId, 

(5.2) 



dvj dvk dvj dvk Ovq dvk 



dtk dtj dtk dtj dtk dt 

which follows from the Lax-Sato equations. We prove this below following the lines 
of Ref. [0. To prove the first symmetry relation, we notice from ( |4.4| ) that dtjVk is the 
constant (z^) part of the Laurent series z^^^dzHj in i.e., the residue Tes{z^dHj). Then, 
using the well known property of residues, we find that this is equal to the constant part 



(w") of the Laurent series z^wdwHj in w. Then, using the equation ( [4.7|) we get: 

(J f) r) nih 

{{z^).w—{z^)^), = -{{zrUw—{z'y)o = {{z^)-W—{z%)o = ^ 

OW OW OW Otj 

This chain of equalities is due to the identity {f-wdwg+)o = —{g+'wdwf-)o = {g-'wdwf+)Q 
for residues of Laurent series. The proof of the other two symmetry relations is similar. 

In fact, the symmetry relations follow from the unitarity condition (|2.2| ) for the 
Schwarz function. To show this, let us outline another proof which does not rely on 
existence of the conformal map. The proof uses only the fact that the variation of the 
Q{z) with respect to the harmonic moments tj, if z belongs to the curve, is purely imag- 
inary (see the end of Sec. 3). Using the expansion of fl in the series ( |3.7| ) and the second 
equality in ( |4.2|) , one concludes that 

/ Hj{z)dHk{z) = (5.3) 
Jcurve 

for all j,k > 0. Therefore, we have: 

dt.Vk = — / z^dHj{z) 
27X1 J curve 



ZmJcuive ZmJcuive \ / I I 



o ■ r I 7;\'"o + y^^tfc^i ^\ dz^ = dt^Vj 
27rz/curve \ 2 ^ 



and similarly for the other derivatives. 

The r-function determines the Hamiltonians as functions of z. Eq. ( [4.4| ) reads 

rr _ , (id' , ^ Z-'^ \ 

' ^ \2dtdt, k dtjdtk] 



The second equation in ( p. 61) gives the formula for the inverse conformal map w{z) through 
the r-function (see (|1.6| ) in the Introduction). Formulas for the coefficients Uj of the direct 
map z{'w) through logr are also available but they have more complicated structure. The 
first two are simple: r = | exp(9j^logr), = d'^logr/dtdti. 
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The r-function itself obeys the dispersionless hmit of the Hirota equation (a leading 
term of the differential Fay identity 0, 0], 0): 

(-0 expf E r;r--X--l=-xpf-E^^-'^Vcexpf-E^H (5.4) 

\n,Tn>l ) \ k>\ ) \ fe>l / 

This is an infinite set of relations between the second derivatives Vnm = c^t™ log^, "i^om = 
9t9f^logr of the r-function. The equations appear while expanding the both sides of ( |5.4| ) 
in powers of z and In particular, the leading terms as C ^ cxd yield 

- - -01 - E = - exp [- E ""-^zA (5.5) 

fc>l V / 

This is in fact the relation H\ = re^° +uo/2 between the Hamiltonians Hi and Hq = logw 
which follows from the definition ( |4.7| ). Substituting ( p.5| ) back to ( |5.4| ), one can eliminate 
f ofc. The resulting relations between Vmn with m, n > 1 is exactly the dispersionless limit 
of the Hirota equations for the KP hierarchy from |H, H], 0). Supplemented by proper 
initial data, satisfying the unitarity condition, equations ( |5.4|) formally solve the local 
problem of moments (or the 2D potential problem), i.e., reconstruction of the Laurent 
expansion of the S^~\ 



The proof of (|5.4|) is similar to the one known in the context of the KP hierarchy 
EE- Let 

oo 

rw = z + "^PjZ'^ (5.6) 

j=0 

be the Laurent expansion of the inverse conformal map w{z), where pj are some coeffi- 
cients. Multiply both sides by C~'^z^~^(w), extract the polynomial part in w and sum 



over k > 1. After some rearrangement and using (5.6) for w(Q one gets the relation 



oo 

r{w{0 - w{z)) E C'-\z\w))^ = rw{z) E C'-\z\w))o 

k=l k=0 

Using ( |4.7| ) and the identity dlogw{z)/dz = Y^kLo z~^~^{z^{w))q, we rewrite it in the 
form 

^1 - wiz) 



Now, plugging the equivalent definition ( [4.4|) of Hn and integrating this equahty with 
respect to (, one finally arrives at the following important relation: 

E^---log^^^(^l^ + logr (5.7) 



Because of ( |1.6D this is the same as 

Let us note, in passing, that the Schwarzian derivative of the inverse map w{z), 

w'{z) 2 \ w'{z) 
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(the classical stress-energy tensor) admits a remarkably simple representation linear in 
logr. Indeed, differentiating both sides of (|5.7|) with respect to z and (, we get 



A simple calculation shows that the limit of the right hand side as C — > 2; is equal to | of 
the Schwarzian derivative (15.81). Therefore, we obtain: 



rw = 6.-i:--'-"|!l^ (5.9) 

k,n>l OtkOtn 

The Hirota equation itself is thus a "splitted" {( 7^ z) analogue of (p.9|) 



6 Conformal maps as a reduction of the dispersion- 
less Toda lattice hierarchy 

The reader familiar with integrable hierarchies of non-linear differential equations is able 
to identify the dynamical system for conformal maps ( [4.5|) , ( [4.6|) with the dispersionless 
limit of the Toda lattice hierarchy 0, 0]. The latter is related with the Whitham hierarchy 
- the theory of solitons with distinct fast and slow variables. The Whitham hierarchy 
appears after averaging over fast variables (see 0, [11| and references therein). The 
dispersionless limit emerges as the genus-zero Whitham hierarchy. Formally, it is a 
semiclassical limit h —>■ of pseudo-differential (or difference) operators. For the Lax 
operator of the 2D Toda lattice, 

00 

L = r(t)e^H + J2 M, (t)e-^'^H (6.1) 

j=0 

one should replace the difference operator e^H by the canonical variable w with the 
Poisson bracket {\ogw,t} = 1. The Lax operator then becomes the Lax function L[w) 
given by a formal Laurent series in w. The Lax function is identified with the conformal 
map z{w) (11.51). The derivatives of the z{w) with respect to the times are given 
by (|4.5| ) which is nothing else than the dispersionless limit of the Lax-Sato equation 
hd/dtjL = [Hj, L], where Hj = {L^)_^_ + ^{L^)q. (The coefficient 1/2 is due to a particular 
choice of gauge in ( |6.1D , where the coefficient in front of the first term is not fixed to be 
1.) The mathematical theory of the dispersionless hierarchies constrained by the string 
equations has been developed in Refs. 0, Q and extended to the Toda hierarchy in 
Ref. 0. For a comprehensive review see, e.g., ]n| . 

In the Toda theory, there are two sets of independent times, tj and ti, and two sets 
of potentials, Uj and Uj. The identification with conformal maps requires them to be 
complex Uj = Uj and the function r to be real. Under the reality 

conditions the semiclassical limit of the second Lax operator, 

00 

L = r{t-h)e-''^t + ^ u,{t)e^^^t (6.2) 
i=o 
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is identified with z{w ^). The reahty condition is consistent with the 2D Toda hierarchy 
and selects a class of solutions relevant to conformal maps. 

This class is reduced to a unique solution by imposing the string equation {z, z} = 1, 
which in the dispersionful case would be [L, L] = h. To clarify the origin of the string 
equation, one needs two more operators. These are the Orlov-Shulman operators 

oo oo 

M = Y^ ktkL'' + t + E ^kL-^ (6-3) 

fc=l k=l 

oo oo 

M = E kikL" + t + E ^kL-" (6.4) 

k=l k=l 

obeying the conditions [L, M] = hL ; [L, M] = —hL. Then the string equation follows 
from the relations L = L~^M, M = M. From Sees. 2,3 it follows that the disper- 
sionless limit of the operator L^^M enjoys a simple geometric interpretation: it is the 
Schwarz function S{z) of an analytic curve. 



7 More connections and equivalences 

Analytic curves and dispersionless hierarchies. Thus one particular solution of the dis- 
persionless 2D Toda hierarchy describes evolution of the univalent conformal map of a 
domain bounded by a simple analytic curve to the exterior of the unit circle. The set of 
times of the hierarchy appears to be equivalent to the set of harmonic moments of the 
domain, whereas the conformal map itself is the dispersionless limit of the Lax operator. 

This proposition can be reversed: univalent conformal maps of the exterior of the unit 
circle generate a solution of the dispersionless Toda hierarchy. The solution is selected 
by the string equation (|3.1|) . 

Other solutions of the dispersionless Toda hierarchy are selected by more general 
string equations. The latter are characterized by any two functions / and g forming 
a canonical pair: {f,g} = f- Then the general string equations consistent with the 
hierarchy are L = f~^{L,M), M = g{L,M) (see 0). We expect that some of them are 
also relevant to the conformal maps of simply connected domains. It is likely that other 
types of string equations describe mappings to or from domains other then the exterior 
of the unit circle and also nonunivalent maps. 

Moreover, we expect that other integrable hierarchies of the hydrodynamic type, for 
instance KP hierarchy and its reductions, also describe certain classes of conformal maps 
and curves others than analytic. We plan to address this question elsewhere. 

The T-function for curves. All this suggests to introduce a general notion of the 
r-function for curves. For simple analytic curves, this is a universal function which 
determines the curve by means of the Laurent series 

oo Q 

x^ + y^-t = J2 T^eiktk {x + iyf + {x + iyy^—logr) (7.1) 
k=i c'tfc 

The dispersionless Hirota equation ( |5.4|) provides a set of differential equations for the r- 
function. Again, it has many solutions. The proper solution is selected by the unitarity 
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condition ( p.T]) and initial data. For instance, if the curve can be reached by small 
deformations of a circle, one sets dt,.logT = at all = 0. 

Moving curves and the c = 1 string theory. There is another intriguing relation 
(in fact equivalence) between deformations of analytic curves and the genus-0 topological 
sector of 2D gravity coupled to c = 1 matter . This follows from the known equivalence 
between the latter and the dispersionless Toda hierarchy restricted by the so-called W^i+oo- 
constraints 0, H, The interpretation of the objects of the genus-0 string theory in 



terms of analytic curves is straightforward. The positive (negative) momentum tachyon 
one-point functions < Tn > are moments of the interior of the domain f„ (vn), the 
partition function of the genus-0 string is the r-function of analytic curves, the Schwarz 
function S{z) is the superpotential and, finally, the l^i+oo-constraints are essentially 
equivalent to the string equation ( ^.1| ) or the relation ( 4.11 ). In fact, they are nothing 
else than the unitarity condition (|2.71). 



In its turn, the c = 1 genus-0 string theory, as well as the relevant solution to the 
dispersionless hierarchy, is known to be equivalent to the planar limit of the two-matrix 
model (see and references therein). This gives a representation of the r-function of 
analytic curves as iV — > oo limit of the NxN two-matrix integral t = J DMDM exptrW , 
with the potential W = EkihM^ + hM^) + MM. We plan to discuss this subject 
elsewhere. 

The interpretation of the genus-0 string theory as a simply minded classical theory of 
potential may turned out to be fruitful for both subjects. It is likely that the 2D gravity 
coupled to c < 1 matter, being described by various reductions of the dispersionless KP 
hierarchy, also enjoys a geometric interpretation. 

Laplacian growth problem. This paper has been stimulated by our interest in the 
Laplacian growth problem (see fl^ for a review). This problem (a source of interest- 



ing mathematics and a great deal of important application) may give further insights in 
the matters we discuss in the paper. It seems to be the place for a short introduction. 
This is the problem of a moving interface between two incompressible liquids with dif- 
ferent viscosities on the plane. Let, say an exterior of a simply connected domain be 
occupied by a viscous fluid (oil), while less viscous liquid (water) occupies the interior. 
Water is supplied by a source at 2; = 0, while oil sinks to a source at the infinity, so the 
interface moves. Experiments and numerical simulations suggest that any smooth inter- 
face, regardless of its initial shape, develops a finger-like pattern with a universal fractal 
characteristics. The hydrodynamics of an ideal interface (with zero surface tension) is 
described by the Darcy law 

I (-) 

where Vn is the velocity of the interface and dp/ dn is the normal derivative of the pressure 
on the interface. The pressure is constant in the water domain and on the interface, while 
in the oil domain it obeys the Laplace equation V^p = with the asymptotic behavior 
at the infinity p — — l/21og|2;|. The latter indicates a sink at the infinity. The Darcy 
law implies [jl8| that all harmonic moments tk of the oil domain are not changed while 
the interface moves, but the area of the water domain grows linearly in time and thus 
can be identified with the time. (For connections with the inverse potential problem see 
||T9| .) The problem then becomes: find the evolution of the domain as a function of the 
area t at fixed moments tk. This evolution is described by the string equation ( |3.1| ). 
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This equation has a long history. It appeared in 1945 in Ref. I^Oj or even earher in the 
mathematical theory of oil hydrodynamics. In our approach, this equation is the basis 
of the symplectic structure of the conformal maps. Applications of integrable hierarhies 
to the Laplacian growth problem are addressed in Ref. |l^ 

At the time of completion of the manuscript E.Ferapontov informed us that J. Gib- 
bons and S.Tsarev discussed a relation between Benney equations and conformal maps 
of slit domains |2ll]. We would like to thank J. Gibbons and S.Tsarev for informing us 



about their recent paper 22 
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Appendix: Ellipse growing from a circle 

Here we demonstrate how the Lax-Sato equations describe an ellipse growing from a 
circle. It is the simplest but nevertheless instructive example which in particular allows 
one to compute the r-function at all tfc = but t, ti, ti, ^2, h- Consider an ellipse with 
half-axes a, b centered at Zq = Xo + iyo and rotated by the angle a: 

(cosa(x - Xq) - sin a{y - yo)^ (sin a{x - xp) + cos a{y - yo)f _ 
a? ~ 

The Schwarz function of the ellipse is 

^2 I ip, n 2ia 

S{z) = e^^" {z - z,) + -Z0- -^^^(z-zoy-e^^'-ia'-b^) 

The Laurent series of the Schwarz function ( p.5|) 

S{z) = 2t2Z + h + - + ^ + ^ + ... 

z z^ z^ 
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gives the moments of the exterior and the interior. The only nonzero moments of the 
exterior are ti, t2 and their complex conjugate: 

2t2 = e^*"^^, h = zo- 2t2Zo, t = ah 
a + 

Contrary, none of moments of the interior vanish. The first two are 

t(ti + 2tit2) t{ii + 2tit2)^ 2tH2 



4t2t2 (l-4t2t2)2 l-4t2t2 



Adding dfVQ = 2 logr, one may check the symmetry relations ( |5.2| ) and find the r-function 
for the ellipse: 

bgr = -tHogt - -t^ - -tHog 1 - 4t2t2 + t — ^ - ^ 

^ 2^4 2^^ ' 1- 4t2t2 



Let us note that this function (for ti = and t2 real) was obtained p3| as the free energy 
of a classical Coulomb gas in an external field. 

The Laurent series for the conformal map from the exterior of the unit circle to the 
exterior of the ellipse is truncated: 



z{w) = rw + Uq + UiW 



The coefficients of the conformal map are: 



4 ' ' 1 - At2t2 4 ' 'I- At2t2 

h + 2tit2 
"° = '° = T^1S 

The first two Hamiltonians are: Hi = rw + ^Uq , H2 = r'^w'^ + 2ruoW + rui + ^Mq. Higher 
fiows deform the ellipse. The Lax-Sato equations plus the string equation and their 
conjugate describe how the ellipse grows from the circle. 
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